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Abstract 

The allowed parameter space for the lightest neutralino as the dark matter is 
explored using the Minimal Supersymmetric Standard Model as the low-energy ef- 
fective theory without further theoretical constraints such as GUT. Selecting values 
of the parameters which are in agreement with present experimental limits and 
applying the additional requirement that the lightest neutralino be in a cosmolog- 
ically interesting range, we give limits on the neutralino mass and composition. A 
similar analysis is also performed implementing the grand unification constraints. 
The elastic scattering cross section of the selected neutralinos on ^^Al and on other 
materials for dark matter experiments is discussed. 



*Subniitted to Astroparticle Physics, 19 Feb. 96 
^Corresponding author, gabutti@vms.mppmu.mpg.de 



1 



1 Introduction 



Particle candidates for dark matter are classified as hot or cold depending on whether 
they were relativistic or not at the time they decoupled from thermal equilibrium. Any 
particle candidate for cold dark matter implies an extension of the Standard Model. At 
present the most interesting candidate is the lightest neutralino of the Minimal Super- 
symmetric Standard Model (MSSM). 

Any weakly interacting massive particle (WIMP) considered as a dark matter candi- 
date is subject to at least two constraints: 1) its relic abundance must be cosmologically 
interesting, say 0.025 < flh'^ < 1, in units of the critical density and 2) its existence must 
be in accord with present experimental limits, provided mainly by the LEP experiments. 

In the standard approach, based in solving the kinetic Boltzman equation, the relic 
abundance of WIMPs is given roughly by : 

nh^ ~ ^ (1) 

where {(Tannv) is the thermal average at freezeout of the annihilation cross section times the 
relative velocity in units of c. This is a very remarkable result: Qh"^ ~ 1 for typical weak 
cross sections. Suppose the annihilation of a WIMP of mass m proceeds via the exchange 
of a particle of mass M (where = xMf) coupled with the strength g [g"^ = ygzuu)^ 
where we have introduced the scaling factors x and y to express the mass M and the 
coupling g in terms of the Z boson mass and its coupling to neutrinos. Using Eq. ([^) and 
the expressions for {(Tannv) which follow from dimensional arguments we have: 



- M ife - S (t^)^ 0.4 X 10-' c.^ f„. » « M (2, 



{(^annv) ~ 1/' ^ ^ = ( ^ X 0.2 X lO^^^ cm^ for m > M 



IGvr^ \lTeV 
and requiring Qh'^ < 1 leads to a generalized Lee- Weinberg bound of: 

- X on GeV) <m<yxO(l TeV) . (3) 

y 

This is a useful qualitative constraint on the mass of a WIMP when used together with 
experimental limits on the parameters x and y (of course x and y are also constrained 
by the requirement that flh"^ should not be too small). It is convenient to distinguish 
two cases: a neutrino-like WIMP which interacts only via Z-boson exchange and a non 
neutrino-like WIMP which can annihilate via exchange of new particles. 

In the first case, the presently available experimental constraint on the "invisible" 
width of the Z, ATinv/Ti, < 0.05 0, gives y < 0.05 (for of course x=l and m < Mz/2). 
Therefore, from Eq. (|^) we obtain m > (9(20 GeV). This is an order of magnitude 
estimate; in particular we neglect the difference between Majorana and Dirac particles. 

The second case is realized in the Minimal Supersymmetric Standard Model with the 
lightest neutralino as the dark matter candidate. Here the analysis is much more involved 
because of the complexity of the model and its large parameter space. A considerable 
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amount of work has been devoted to the neutrahno as a dark matter candidate How- 
ever in most cases several additional theoretical assumptions are used (such as radiative 
electroweak symmetry breaking and universal boundary conditions at the GUT scale for 
the parameters of the soft breaking of supersymmetry) which may be too restrictive and 
certainly go beyond the MSSM as the low-energy effective theory. 

In view of the dark matter search experiments, we explore in Sec. ^ the most general 
scenario for neutralinos as dark matter within the MSSM as the low-energy effective 
theory without any further theoretical constraints such as the grand unification constraints 
(GUT). For completeness, the obtained neutralino masses and compositions are compared 
with the results derived using the GUT constraints. 

In Sec. ^ we consider the direct detection of neutralino dark matter via elastic scat- 
tering on nuclei. The dependence of the axial coupling (spin dependent) on the nuclear 
model and on the quark spin contents is discussed for selected neutralino compositions. 
In Sec. § we calculate the neutralino cross section on some of the nuclei presently used or 
planned to be used for dark matter detection. We analyze the low mass region and give 
the dependence of the cross section on the neutralino mass and composition. Prospects 
for the detection of light neutralinos are discussed in Sec. |^. 

Our analysis differs from previous work on neutralino detection ||^ because our more 
general assumptions allow the possibility of having low mass neutralinos as dark matter 
candidates and we evaluate their cross section on nuclei. In the present work neutralinos 
with relic abundance below the cosmological bound Qh!^ < 0.025 are not considered as 
dark matter candidates. This is different from the approach used in other work [^] where 
neutralinos with relic abundance below the cosmological bound are still considered as 
dark matter candidates and their contribution to the dark matter halo is evaluated by 
rescaling the local dark matter density. 

2 Neutralino in the MSSM as dark matter candidate 

In the present section we address the question of the most general limits on the lightest 
neutralino mass and its composition which follow only from the two constraints 0.025 < 
fih"^ < 1 and consistency with the present data from accelerator experiments, without 
any further theoretical assumptions. As we shall see, with this approach it is possible to 
obtain quite strong qualitative conclusions. 

The stable neutralino is the lowest mass superposition of neutral gauginos and higgsi- 
nos: 

Xi = ZuB + ZuW' + Z.^Hl + Z,,Hl (4) 
or, in the photino, zino, higgsino basis: 

Xi = a^ + hZ + Z^^Hl + Zi4^2° (5) 

with 7 = cosOwB + smOwW^ , Z = — smOy/B + cosOwW^ and 6w the Weinberg angle. 
The neutralino composition is defined by the neutralino mass matrix, which in the basis 
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of Eq. (D is: 



/ Ml 


—Mz cos (3 sin 9w 
V Mz sin P sin 6*1^ 





Ms 

Mz cos /3 cos 6*11/ 
—Mz sin /3 cos6'vv 



— Af^ cos/? sin 6' 
M^ cos /5 cos 9 




V M^ sin P sin 6'vi/ \ 
'w — Mz sin /3 cos 6'i4/ 
-/i 

y 



(6) 



2.1 Scan of the parameter space 

The neutralino mass matrix depends on several, in general free, parameters of the 
model: the gaugino masses Mi and M2, the higgsino mass parameter /i, and tan j3=V2/vi, 
where vi and V2 are the vacuum expectation values of the two Higgses present in the 
model. All of these are independent parameters of the low-energy effective lagrangian. 
In the present study of the unconstrained low-energy effective theory, the values of these 
parameters are chosen randomly in the ranges listed in Tab. |l|. The scanning procedure is 
such that log(M2), log(/i), log(tan/3) and the ratio M1/M2 are chosen with flat probability. 
We restrict our analysis to neutralinos lighter than 100 GeV. 

Calculations with the grand unification constraints are also performed using the same 
ranges but with Mi = |M2 tan^ ^vi/ ~ O.5M2 and with M2 ~0.3mg. We use for the gluino 
mass rrig the experimental lower bounds from CDF and DO ^ resulting in M2 >50 
GeV. It should be noted that the present experimental limits on the gluino mass are 
derived only for a specific choice of the MSSM parameters. 

We scan the parameter space incorporating the following existing experimental limits 
from accelerator experiments. 1) The limit on the Z boson invisible width is T^nv < 
8.4 MeV 10]. 2) Heavier neutralinos are unstable and could be observed in e~^e~ collisions 
via their decay products for the range of x masses which are kinematically accessible. We 
use the mass-dependent limits on the coupling constants given in Ref. [0. 3) The limit 
on the chargino mass is M^± > 45 GeV . The chargino mass matrix is a function of the 
same parameters: 

M2 V2Mw sin p 

V2Mw cos P fj, 

and does not depend on the ratio M1/M2. 

The scanned parameter space is shown in Fig. |l| with and without grand unification 
constraints after incorporating the limits from accelerator experiments. The empty regions 
in the M2 — /i plane are excluded by the experimental limits. The quoted chargino mass 
limit is at present by far the strongest constraint on the parameter space as shown in Fig. |^ 
where the contour line M^± = 45 GeV is plotted for tan/3=l and tan/5=50. Without 
grand unification constraints, the region corresponding to small values of M2 and \fj,\ is 
excluded for tan /5 > 3. 

2.2 Relic abundance 

We then calculate the relic neutralino abundance in the allowed parameter region. We 
solve the Boltzman equation following Ref. [§, using a complete set of supersymmetric 
annihilation amplitudes 0. For this we need the squark, slepton, and Higgs masses. We 
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assume all squarks have the same mass Mgq and all sleptons the same mass Mgi- Since 
the purpose of our study is to explore the parameter space, we need limits which are 
model-independent, whereas most experimental limits are based on certain assumptions 
on the values of the parameters. Thus we have chosen to repeat our calculation for several 
different values of the squark, slepton and CP-odd Higgs masses as shown in Tab. |^. In 
the case of sets A-C, the upper bound in the neutralino mass is given by the slepton mass 
since we require < Mgi. 

The lowest values of the squark masses are taken as 45 GeV, since lighter particles 
are ruled out by LEP data on the Z. CDF has shown that the squark mass limits 
obtained from pp collisions are dependent on the assumed squark decay modes and that 
for many configurations of the MSSM parameters the squark mass is unconstrained. A 
recent report by DO P] gives limits on the squark mass as function of the gluino mass for a 
particular set of MSSM parameters, assuming that all squarks have the same mass. Even 
for this particular set of parameters, squarks can be as light as ~50 GeV for rrig >300 GeV. 
OPAL [|TU| has set limits of about 45 GeV specifically for the t mass, which may be 
significantly lighter than the other squarks due to the left-right mixing. 

For sets A-C the CP-odd Higgs mass Ma is randomly chosen in the range 25-70 GeV 



in the regions of the plane Ma-Mh allowed by present LEP data [11|. In sets D and E 



Ma is fixed at 200 GeV and 1 TeV respectively. The other Higgs masses are calculated 
from Ma and tan (3 using radiative corrections [|12| which introduce a dependence on the 
top quark mass Mt and Msq. The top quark mass is set to 170 GeV. 

The result of our calculation for the relic abundance Vth'^ obtained without the grand 
unification constraints are plotted in Fig. |^ for sets B-E. The cosmologically relevant 
neutralinos are selected imposing the condition 0.025 < VlK^ < 1. It must be noted 
that we do not rescale the local dark matter density for neutralinos with relic abundance 
below the cosmological bound Vth"^ < 0.025 since we do not consider them as dark matter 
candidates. 

For low squark and slepton masses (set B) the condition VLh? < 1 is obtained for 
neutralino masses above 2 GeV. Decreasing the squark masses to 45 GeV (set A) gives a 
relic abundance very similar to the one of set B and the same value of 2 GeV for the lowest 
neutralinos, which can then be regarded as an absolute lower bound on the neutralino 
mass in the MSSM. This is one of the main conclusions of this study. The lower bound 
in the neutralino mass rises with Mgq and Mgi to reach ~20 GeV in sets D and E. 

The relic abundance evaluated using the grand unification constraints is plotted in 
Fig. ^ for sets B-E. Here neutralinos have a lower bound on the mass of 20 GeV, due to 
the experimental constraints on the gluino mass. A relic abundance VLh"^ << 1 is obtained 
for sets A (not shown) and B where we use low squark and slepton masses. 



2.3 Neutralino composition 

In Fig. ^ we show the gaugino purity + separately the bino purity 

for the cosmologically relevant neutralinos derived without grand unification constraints 
combining data sets A-E. The points at bino purity Zl^ = cos^ 9w — 0.77 correspond 
to pure photinos. It is seen that the cosmologically relevant neutralinos are dominantly 
gauginos (mostly bino and photino), with a substantial higgsino component coming in 
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above 30 GeV. The neutralino composition calculated with GUT is shown in Fig. ^ and it 
is similar to the composition obtained without grand unification constraints except that 
the photino population of the gaugino-like neutralinos is strongly reduced. The higgsino 
composition is shown in Fig. ^ without grand unification constraints. Since we consider 
values of tan/? >1 we have Zf^ > Zf^ with Z13 = —Z^^ when tan/3=l. The higgsino 
composition with grand unification constraints is similar to the one shown in Fig. ^ 

The depopulation in the region with large gaugino- higgsino mixing is due to the scan- 
ning of the allowed parameter space and to the requirement of having cosmologically 
relevant neutralinos. The neutralino is a mixed state only when the smaller of the two 
gaugino masses, M2 or Mi, has a magnitude similar to This is a fine-tuning which is 
not included in our scanning procedure, in consequence we have a substantial reduction 
in the population of mixed neutralinos. Especially for large M2 there are few points in the 
parameter space corresponding to mixed neutralinos. In this case, for values of Mi not 
much smaller than M2, mixing is allowed for large values of giving neutralinos heavier 
than our upper bound in the neutralino mass of 100 GeV. The only mixing allowed for 
large M2 corresponds to both Mi and being small. 

The dependence of the relic abundance on the neutralino composition is shown in 
Fig. 1^ for sets B and D without grand unification constraints. Pure gauginos do not have 
couplings to the Z boson and have higher relic abundance than pure higgsinos. The relic 
abundance of pure photinos is bigger than that of pure bino because the former does 
not couple to the Higgs bosons. At large tan/5 the annihilation is dominated by Higgs 
exchange which is proportional to the gaugino-higgsino mixing and goes to zero for pure 
states. As a result, large gaugino-higgsino mixtures have low relic abundance and are 
excluded by the request of having > 0.025. The dependence of the cosmologically 
relevant neutralino compositions on tan /3 is shown in Fig. ^ for sets B and D. Gaugino- 
higgsino mixtures are relevant only at small tan /? while at large tan /5 our cosmologically 
relevant neutralinos are mainly pure gauginos with exception of set E which has a very 
heavy CP-odd Higgs boson. There is a substantial photino component for tan (5 <2 while 
at large tan/3 the photino component is strongly reduced due to the structure of the 
neutralino mass matrix. Due to the large squark mass, photino have relic abundance 
f2/z^ > 1 in the case of sets D and E. 

3 Neutralino detection 

Dark matter neutralinos can be directly detected via their elastic scattering on nuclei. 
The scattering cross section has two components. An effective axial-vector interaction 
gives a spin- dependent (SD) cross section which is non-zero only for nuclei with net spin. 
Scalar and vector interactions give spin-independent (SI) cross sections which involve the 
squares of nuclear neutron and proton numbers. The relative strengths of the two parts 
of the neutralino interaction on nuclei depends on the neutralino composition. A pure 
gaugino couples only to squarks and sleptons and has only a SD interaction. The sleptons 
are not relevant for interactions with nuclei. A pure higgsino couples only to the axial- 
vector part of the Z boson and thus has only a SD interaction. (The Yukawa coupling 
to squarks is negligible for scattering on the proton.) Mixed states have both SD and SI 
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interactions, the later is due to squark exchange and Higgs exchange and is proportional 
to the zino-higgsino mixture. 

The actual scattering cross section depends on the coefficients Z\j that specify the 
neutralino composition, and on the masses M^^ and of the exchanged particles. To 
obtain the neutralino scattering on nuclei for each set of supersymmetric parameters, we 
first calculate the coefficients of the effective 4-fermion operators describing neutralino- 
quark scattering. These operators are used to give the neutralino-nucleon amplitudes, 
which are then related to the nuclear matrix elements. 

We shall focus on calculating the cross section for neutralino scattering on a nucleus 
at zero momentum transfer. This incorporates the essentials of the physics of neutralino 
cross section on nuclei. In Sec. ^ we will include in the evaluation of the recoil energy 
spectrum and detection rate nuclear form factor effects due to the momentum transfer. 
In the following, we will call the integral of the zero momentum transfer cross section over 
all momentum transfer Gq. The cross section is derived from: 



(To = 4^ 5a(0) + Ss{^) 



TT 



where the effective axial- vector current 5*^(0) and the scalar current 6*5(0) coefficients 
are calculated at zero momentum transfer (5'=0) and M^ed is the reduced mass. The 



coefficients Sa and Ss are evaluated using the expressions derived in Ref. |]T3 



3.1 Spin Dependent cross section 

The evaluation of the SD cross section is problematic because it depends on the model 
used to describe the spin structure function of the nucleon and on the nuclear model used 
to derive the total proton and neutron spins. In this section we discuss the resulting 
uncertainty in the determination of the axial- vector current coefficient: 

5^(0) = 8J(J + + ^IL^^I' (9) 

with: 

S(n)= [E^.^QJ,(,) (10) 

where J is the nuclear spin and < S^(n) > are the total proton and neutron spins. The 
terms Ag refer to the coupling of the up, down and strange quarks weighted with the 
nucleon quark spin content coefficients AQg. The entire dependence on the nature of 
the neutralino, and therefore of our cosmological considerations, enters only through the 
coefficients Ag. For the determination of < S'p(„) > and AQg, we rely entirely on previous 



work [|lj]- 



It is illustrative to consider pure photino and pure bino interactions for which the 



proton axial coupling coefficients are |13] 



tan^ e^^^AQ^ + ^(AQ, + Ag,)| pure bino (11) 



M2 

sq 



2^sin'^H/|^Ag„ + ^{AQd + Ag,)| pure photino 

sq 
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Since the quark spin coefficients AQg are evaluated for scattering on a proton, the neutron 
axial coupling coefficient is obtained by interchanging the values of AQu and AQd 
according to the different quark structure of the two nucleons. 

In Tab. ^ we have calculated ap(„) from Eq. ( ]TT| ) using for AQ^ the values derived 
from the Naive Quark Model (NQM) the European Muon Collaboration (EMC) 
measurements and from two analyses of the present data on polarized lepton-nucleon 
scattering Global Fit-1 from Ref. |]l6l and the more recent Global Fit-2 from Ref. ||17 . 



The axial coupling coefficients Op are roughly a factor of two larger than a„ due to the 
different quark structure of the two nucleons which enhances the axial coupling on protons 
and therefore on proton-odd nuclei. 

The dependence of the axial-vector current coefficient Sa{0) on the model used to 
describe the nuclear structure is shown in Tab. § where we have calculated the axial- 
vector current coefficient for pure photino and pure bino interactions on the proton-odd 
nucleus ^^Al and on the neutron-odd nucleus ^^Ge. We compare the predictions from the 
Odd Group Model (OGM) |T§ and the detailed Shell Model reported in Ref. for ^^Al 
and in Ref. for ''^Ge. We also show the effects of different model for the spin quark 
coefficients AQg. 

If the OGM is used to describe the nuclear structure, only the unpaired nucleon 
contributes to the axial coupling. In this case the axial- vector current coefficients ^^(O) 
are proportional to for proton-odd nuclei and for neutron-odd nuclei. In the OGM, 
where only the valence nucleon plays a role, the cross sections on nuclei obtained with 
the NQM are roughly a factor of two larger on proton-odd nuclei than those for other 
quark models directly reflecting the cross sections on single nucleons. On the contrary, 
the SD cross section on neutron-odd nuclei derived with the NQM is more than an order 
of magnitude smaller than the one obtained using the other AQg values due to the small 
value of a„ in the NQM. 

It should be noted that the values for the total proton and neutron spins derived in 
the detailed Shell Model calculations for ^^Al have an uncertainty of roughly 30% related 
to the quenching of the spin matrix elements. Quenching is more important in ^^Ge 
and we use the quenching factor Q = 0.833 in the isovector piece of the axial coupling 
coefficients ap(„) as suggested in Ref. ||20[. In this case the axial coupling coefficients have 



to be replaced with a'^ = O.QlTop + 0.0835a„ and a'^ = 0.0835ap + 0.917a„ where ap(„) are 
evaluated in the unquenched case (Eq. (|1T])). The cross section on ^''Al evaluated using 
the detailed Shell Model calculations are roughly 2 times larger than the results obtained 
with the OGM for all the four model used to describe the spin structure of the nucleon. 
The detailed Shell Model gives almost 4 times larger SD cross sections on ''^Ge than the 
OGM for all the three quark models with non zero strange quark content. The axial 
current coefficients for pure photinos and binos are smaller for interactions on '''^Ge than 
on ^^Al mainly because of the different quark structure of neutrons and protons. 

In general neutralinos are mixtures of gauginos and higgsinos. For pure higgsinos, the 
axial coupling coefficient Qp derived from Ref. [|l^ is: 

J (Zl^ - Zl) [Ag„ - AQ, - AQ J (12) 
Z^M^AQ^, + tan^ /?[MjAQd + M^AQ^ 
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1 1 1 
'2sin2/3M2^ 
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where Mu^d,s are the quark masses. The neutron axial couphng coefficient a„ is obtained 
by interchanging the values of AQu and AQd- For the same higgsino composition, the 
difference between Op and a„ is mainly due to the first term of Eq. ([I2|) which is always 
positive for protons and negative for neutrons since Zf^ > Zf^ (see Fig. ^. The second 
term of Eq. (|12D does not change significantly going from to a„ because it is dominated 
by the contribution of the strange quark which is much heavier than the u and d quarks. 
As a result, \ap\ > |a„| and higgsino interactions have the largest cross section on neutron- 
odd nuclei. 

3.2 Spin Independent cross section 

We evaluate the scalar current coefficient Ss from: 

SsiO) = 2A'[Ps\s + Y.Ph^hf (13) 

c,b,t 

where A is the atomic number. The terms Pi are determined by the couplings in the 
effective XXQQ lagrangian, with contributions from squark exchange and Higgs exchange of 
the lightest and the heaviest Higgs bosons. The supersymmetric content of the interaction 
is represented by the Pi terms. The coefficients Aj are given by the matrix elements of the 
quark operators -^^qiQi taken between the nucleon states and are by standard methods 
expressed in terms of the sigma-term measured in pion-nucleon interactions and the 
nucleon mass. For Xi and P|] we use the values given in Ref. ITSll. 



The contribution of the light u and d quarks to the scalar neutralino-nucleon coupling 
can be neglected and the term P^Aj accounts for the coupling to the strange quark (index 
s) and to the heavy quarks c, b and t (index h). Since we do not consider the contribution 
of the u and d quarks, we do not have to distinguish between scalar cross section on 
neutron or proton. 



It must be noted that the SI cross section evaluated using the PjAj terms of Ref. |l3 
is four times smaller than the SI cross section derived in Ref. 0]. This factor of four is 
independent of the neutralino composition and mass. 



4 Cross section on nuclei 

In this section we derive the SD and SI cross sections on some of the material presently 
used or planned to be used for dark matter detection. In contrast to our illustrative 
examples in the previous section, we will now consider arbitrary neutralino compositions, 
as given by our scan for cosmologically relevant cases. 



4.1 27A1 



We start with the proton-odd nucleus ^^Al (100% natural abundance). Dark matter 
detectors made of sapphire (^4/203) are presently under preparation by the CRESST |21 
and EDELWEISS collaborations m 



^The original notation of Ref. |Tj 



is Sq instead of Pi 
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In order to evaluate the SD cross section we chose to use the quark spin coefficients AQg 
derived in Ref. (Global Fit-2) and the results of the detailed Shell Model calculations 
reported in Ref. |jl9| for the nuclear physics. In Fig. ^ the SI cross section is plotted versus 
the SD part for the cosmologically relevant neutralinos without the grand unification 
constraints. There is a 5 orders of magnitude spread in the values of the SD cross section 
for the neutralino masses and compositions considered in our analysis. The spread is even 
more pronounced in the case of the SI cross section. 

The dependence of the SI cross section on the neutralino composition is shown in 
Fig. [ly for data sets B and D and tan/? <2. Since the scalar coupling is proportional 
to the zino-higgsino mixture, pure neutralino states have zero SI cross section. Mixtures 
with a large photino component have very small SI cross section < 10~^ which is almost 5 
orders of magnitude smaller than the corresponding SD part. The cross section is reduced 
to values below 10^^ pbarn for binos and higgsinos with purity > 0.99. 



The dependence of the SI cross section on tan /? is shown in Fig. 0. Although the SI 



cross section is expected to increase with increasing tan /?, the maximum value of the SI 
cross section for our selected neutralinos changes only slightly with tan (3. This is related 
to the dependence of the neutralino composition on tan/3 as discussed in Sec. 0|. At 



small tan (3 gaugino-higgsino mixtures are cosmologically relevant giving the maximum SI 
cross section. At large tan (3 our selected neutralinos are mainly bino or higgsino with a 
reduced SI cross section. 

The SD and SI cross sections obtained without the grand unification constraints are 
plotted versus the neutralino mass in Fig. 12 for the combined data sets A-E. The bands 
of points in the low mass region of the SD cross section are due to pure photinos and pure 
binos as shown in Fig. |13[ The SD cross section of pure binos is ~2 times smaller than 
the one for pure photinos as shown in Tab. ^. The photino contribution to the SD cross 
section is relevant only for small tan/3 because the photino population in our selected 
neutralinos is strongly reduced at high tan (3. The highest SD cross section is obtained for 
the pure gauginos of set A. In this case squark exchange gives a significant contribution 
because of the small squark mass Msg=45 GeV. 



In Fig. 14 the SD cross section is plotted versus the bino purity for sets B and D 



without grand unification constraints. For heavy squarks (set D) the highest SD cross 
section is obtained for gaugino-higgsino mixtures with a substantial higgsino component. 
This is because squark exchange is suppressed and Z exchange becomes dominant. 

The SD cross section for our neutralinos is sometimes even larger than the cross 
section for a massive Majorana neutrino interacting via Z exchange; this happens when 
the squark mass is low and its exchange gives a significant contribution to the cross section. 
In contrast, the maximum values of the SI cross section are two orders of magnitude below 
that for Dirac neutrinos. 

The SI and SD cross sections evaluated with the grand unification constraints are 
plotted in Fig. ^ for the cosmologically relevant neutralinos of the combined data sets 
A-E. The main contribution to the cross section is given by sets C-E since for sets A 
and B the neutralino relic abundance is VlK^ << 1 (see Fig. The horizontal bands 
in the SD cross section are due to the bino contributions of sets C and D since in the 
case of GUT the photino population is suppressed as shown in Fig. |. Comparing the 
cross sections obtained with and without grand unification constraints, we see that the 
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maximal value of ao and the dependence of the cross section on the neutralino mass are 
almost the same. In the case of GUT, the cutoff in the neutralino mass is at ~20 GeV as 



discussed in Sec. 2.2. 



4.2 23Na and ^^^I 

We turn now to the evaluation of the neutralino cross section on the proton-odd 
nuclei ^^Na and ^^''I both with 100% natural abundance. Nal(Tl) scintillator detectors 
are presently used for dark matter searches by several groups ||23|| . 

The cross sections ctq (SI+SD) on ^^Na and on ^^^I are plotted in Fig. ITB| for the 



combined sets A-E without grand unification constraints. Since ^"^Al and ^^Na have similar 
atomic numbers the SI cross sections in the two materials do not differ significantly. Due 
to the difference in the atomic numbers, the SI cross section on ^^^I is roughly two orders 
of magnitude larger than on ^^Na. 

The dependence of the SD cross section on the neutralino composition is the same for 
Na, I and Al because these nuclei have an unpaired proton. The dependence of the SD 
cross section on the absorber material is given by the reduced mass, the nuclear spin and 
the total proton and neutron spins < 5'p(„) >. In the case of ^^Na and ^^^I detailed Shell 



Model calculations are not available and we use the Odd Group Model where only 
the unpaired proton contributes. The resulting axial-vector current coefficients Sa for 
pure bino and pure photino (purity=l) obtained using the quark spin coefficients Global 
Fit-2 are listed in Tab. ^ The SD cross section on ^^Na is almost three times smaller 
than the values derived for ^''Al with the OGM and the spin quark coefficients Global 
Fit-2 (see Tab. This is due to the different total proton spins and nuclear spins. In 
the case of ^^''I, the SD cross section evaluated with the Global Fit-2 is almost two orders 
of magnitude weaker than the SI part and can be neglected in the evaluation of ctq. 



The bands of points in the total cross section of Fig. [T^ are due to the axial coupling 
of neutralinos with a large photino components, for which the SI cross section is many 
order of magnitude weaker than the SD part. 

4.3 ^^Ge and ^^^Sn 

We will discuss now the case of neutron-odd nuclei starting with ''^Ge (7.8% natural 
abundance). Although the Ge dark matter detectors presently used are made of 
natural germanium which is mainly ^^Ge (36.5%), ^^Ge (27.4%) and ™Ge (20.5%), the 
use of enriched ^^Ge detectors is planned for the future p5[| . 

In order to evaluate the SD cross section we use the quark spin coefficients AQg 
derived in Ref. (Global Fit-2) and the results of the detailed Shell Model calculations 



reported in Ref. |2y] for the nuclear physics. In Fig. |1^ the SD cross section is plotted 
versus the bino purity for data sets B and D without grand unification constraints. In 
both cases the highest SD cross section is obtained for pure higgsinos or for mixtures 
with a substantial higgsino component. In the case of set B this is different from what 
we calculated for Al (see Fig. |l^ where the maximum SD cross section was obtained for 
gaugino-like neutralinos. This is due to the different values of the axial coupling coefficient 



of pure higgsino in the case of proton or neutron odd nuclei as discussed in Sec. 3.1. 
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Since the contribution of the hght u and d quarks is neglected in the computation 
of the SI cross section as discussed in Sec. |3]^, the dependence of the axial coupling on 
the neutralino composition does not differ from the one discussed for the proton-odd ^^Al 
nucleus. In Fig. |T8| the SI cross section is plotted versus the SD part for data sets B and 
D without grand unification constraints. The SI cross section is generally bigger than the 
SD part because the quite large value of the atomic number enhances the scalar coupling. 
For a large selection of neutralino compositions, the SD part can then be neglected in the 
evaluation of the total cross section. 

In Fig. ^ the SD and SI cross sections are plotted versus the neutralino mass for the 
combined sets A-E without the grand unification constraints. The bands of points in the 
low mass region of the SD cross section correspond to the contribution of pure photinos 
of sets A-C with set A giving the highest contribution. The structure due to pure bino 
is less evident than in the case of Al. It must be noted that the photino population 
decreases with increasing values of tan /3 as discussed in Sec. p73| . In Fig. |20| we plot the 
SD and SI cross sections for the combined sets A-E with grand unification constraints. 
The dependence of ctq on the implementation of the grand unification constraints is the 



same as discussed for Al in Sec. 



Another material planned for dark matter detectors is Sn [Q. Natural tin consists of 
different isotopes among which the only odd-isotope with a relevant natural abundance 
is ^^^Sn (7.7%). The axial- vector current coefficients for pure bino and pure photino 
interactions on the neutron-odd isotope ^^^Sn are listed in Tab. |^. The axial coupling of 
pure binos and pure photinos evaluated with the OGM [18| and the Global Fit-2 is roughly 
three times larger on ^^''Sn than on ''^Ge mainly because of the different nuclear spin J. 
The SI part can be derived from the SI cross section evaluated for ^''Al with appropriated 
scaling due to the different nuclear mass. The cross section ctq (SI-t-SD) on ^^^Sn is plotted 
in Fig. pri for the combined sets A-E without grand unification constraints. Due to the 
large value of the atomic number, the SI cross section is generally bigger than the SD part 
and, for a large selection of neutralino compositions, the axial coupling can be neglected 
in the evaluation of ctq- 



5 Prospects for detection of low mass neutralinos 



The present experimental limits for WIMP dark matter are given by Ge [24| and 



Si [E7| ionization detectors and Nal scintillators P3. Due to their relatively high recoil- 



energy threshold, such detectors are not sensitive to low mass (M^ < 10 GeV) dark matter 
particles. 

In order to explore the low mass window, a dark matter search experiment (CRESST) 
based on the use of a sapphire (AI2O3) cryogenic detector with a low energy threshold is 
currently under preparation [^. The first stage of this experiment (CRESST- 1) will use 
a 1kg sapphire detector with an expected energy threshold of 0.5 keV and a full width 
half maximum energy resolution A£'=200 eV. 

We perform a rough estimate of the CRESST-1 sensitivity to neutralino dark matter 
assuming a fiat distribution of the radioactive background at 1 count /keV/kg/day. A 
discussion of different radioactive background sources can be found in Ref. For 
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each dark matter particle mass, the simulated background data is fitted with the sum 
of a flat component for the background and the calculated recoil energy spectrum for 
^^Al. The exclusion limit is defined as the cross section for which the fit gives a 
value corresponding to the 90% confidence level. The expected recoil energy spectrum is 
calculated using an exponential form factor to account for the loss of nuclear coherence at 
high momentum transfer . We use the same exponential form factor for both SD and 
SI interactions. Due to the small nuclear radius the effect of the nuclear form factor on 
the shape of the recoil energy spectrum is negligible for Al. We use a Maxwell-Boltzman 
velocity distribution for the dark matter halo with an average velocity of 270km/s and 
an upper cutoff at the escape velocity of 575km/s in the rest frame of our galaxy. The 
local halo density depends on the model used to describe the structure of our galaxy with 
a resulting uncertainty in the detection rate of roughly a factor of two f^. In this work 
the density of dark matter particles is assumed to be 0.3GeV/cm^. 

In Fig. 1^ the cross section ctq (SI + SD), calculated without the grand unification 
constraints for the combined data sets A-E, is compared to the rough estimate of the 
expected sensitivity of the first stage of the CRESST experiment assuming a measurement 
time of 1 year. The expected low energy threshold of the cryogenic detector used in the 
CRESST experiment combined with the low atomic number of Al provide an appreciable 
sensitivity to dark matter particles starting at ~1 GeV. It is important to note that ^^O 
does not contribute to the SD cross section on sapphire because of the even number of 
nucleons. Due to the small nuclear mass, the SI cross section on ^^O is roughly an order 
of magnitude lower than on ^^Al and can be neglected in the evaluation of the sensitivity 
of a sapphire detector. 

For completeness we have plotted in Fig. |23| the normalized interaction rate in units 
of counts/kg/day/pbarn that is, we show the interaction rate assuming a cross section of 
1 pbarn. The measured rate in a given detector will be less than the interaction rate, due 
to threshold and efficiency effects. The interaction rate is evaluated using the analytical 
expression derived in Ref. with a zero energy threshold^. The variation with mass 
reflects the change in the dark matter flux combined with the effect of the nuclear form 
factor which is relevant for scattering of heavy neutralinos on large nuclei. The variation 
with nucleus reflects the difference in the number of target nuclei per unit weight. In the 
case of ''^Ge and ^^^Sn the interaction rate is calculated for a 100% enriched detector. 
We have separated the contribution of ^^Na and ^^''I to the total interaction rate of 
Nal detectors (50% Na and 50% I atoms). In the case of sapphire detectors we have 
considered only the rate due to neutralino interactions on the ^''Al content of sapphire 
since the neutralino cross section on ^^O is negligible. 

A final question concerns how the counting rate for a given neutralino will vary when 
we change the target material. This is of interest in connection with two points. One is 
the verification of a presumed dark matter signal and the other is discovering the nature 
of the neutralino. Should there be an indication of a dark matter signal in a detector 
it will be of course of the utmost importance to verify that it is a true signal and not 
the result of noise or backgrounds of some kind. An aid in doing this will be variation 
of the target material, leading to changes which are not necessarily those of the noise or 

^We use the expression given for {Ratcc) 
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background. Secondly, varying the type of nucleus can help to determine the character of 
the particle, to see if its interactions are dominantly SI or SD, if it is stronger on neutrons 
or protons and so forth, thus helping to pin down the type of neutralino. Disentangling 
and understanding these various effects will be a question for much further detailed study. 
However, we can give a first impression on the basis of the calculations in this paper by 
presenting examples with some of the various candidate neutralinos that we have found 
in our study. 

In Fig. ^ we show the dependence of the cross section ctq for different neutralino 
compositions on various materials in the case of set B with tan (3=2, Ma=50 GeV and 
the Higgs mixing angle a = —1.24. We use Global Fit-2 for the quark spin coefficients 
and the OGM for the total proton and neutron spins for Na, I and Sn. For Al and Ge 



we use the shell model as discussed in Sec. ^Tl|. The highest cross section is given by 



zino-higgsino mixtures since pure higgsino and gaugino states do not have SI interactions, 
which increase with A^. For other examples, which we take to be pure higgsino or gaugino, 
we show no entries for ''^Ge since it has spin zero and thus there is no SD interaction. 
For pure higgsinos we have higher cross sections on neutron-odd nuclei, as discussed in 
Sec. PTT| . The cross section for pure higgsinos and pure photinos on proton-odd nuclei 
tend to have roughly the same value. In the case of pure higgsino or gaugino states, ctq 
does not scale with since the SI part is absent and the dependence on the materials is 
due to the nuclear model. It should be stressed that Fig. ^ represents only a particular 
case. As discussed in the previous sections, the cross section on nuclei depends not only 
on the neutralino composition but also on the squark and Higgs masses and tan j3. 



6 Summary 

The minimal supersymmetric standard model as the low-energy effective theory has 
been used to explore the most general scenario for neutralino dark matter with the present 
accelerator data and the requirement 0.025< Qh"^ <1 as the only constraints without 
further theoretical assumptions. In order to obtain limits on the neutralino mass which 
are model-independent, we scanned the parameter space randomly choosing the values of 
the free parameters of the model using different sets for the squark, slepton and CP-odd 
Higgs masses. 

We have found an absolute lower bound for the mass of the cosmologically relevant 
neutralinos of > 2 GeV. The cosmologically relevant neutralinos are dominantly gaug- 
inos with an higgsino component coming in for masses above 30 GeV. We performed a 
similar analysis including the grand unification constraints and obtaining dark matter 
neutralinos with a lower bound on the mass of 20 GeV. 

For the cosmologically relevant neutralinos, we have calculated the spin dependent 
(SD) and spin independent (SI) cross sections on different nuclei under consideration for 
dark matter detection. We have found a 5 orders of magnitude spread in the values of the 
SD cross section for the neutralino masses and composition considered in our analysis. 
The spread is even more pronounced in the case of the SI cross section. 

We have discussed the dependence of the SD cross section on the nuclear model and 
on the nucleon quark spin content. We have used the Odd Group Model and the detailed 
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Shell Model calculations to evaluate the total proton and neutron spins and have found 
that the uncertainty resulting from the choice of the nuclear model is within a factor of 
two for the quark spin models with a non zero strange quark contribution. The naive 
Quark Model, where only up and down quarks are considered, gives the smallest cross 
section on neutron-odd nuclei and the highest axial coupling for proton-odd nuclei. 

For our selected neutralinos, the contribution of the photino component to the SD 
cross section is relevant only for small tan (5 because the photino population is strongly 
reduced at large tan /3. The highest SD cross section is obtained for the pure photinos 
and pure binos of set A for which squark exchange gives a significant contribution be- 
cause of the small squark mass Msg=45 GeV. For heavy squarks the highest SD cross 
section is obtained for gaugino-higgsino mixtures with a substantial higgsino component 
because squark exchange is suppressed and Z exchange becomes dominant. The higgsino 
contribution to the SD cross section is more relevant in neutron-odd nuclei. 

In the case of ^^Al, the SD and SI cross sections are comparable in magnitude for 
a large selection of neutralino compositions. The SI cross section on ^^Ge is generally 
bigger than the SD part because the larger value of the nuclear mass enhances the scalar 
coupling. For a large selection of the parameter space, the SD part can be neglected in 
the evaluation of the cross section gq. This is true also for ^^'^I where the SD cross section 
is almost two orders of magnitude smaller than the SI part and can be neglected. Due to 
the different values of the total proton spins and nuclear spins, the SD cross section on 
^^Na is almost three times smaller than on ^^Al. 

The SD cross section increases with the increasing higgsino component due to the Z 
exchange contributions. SI cross sections are large for neutralinos which are zino-higgsino 
mixtures. Large cross sections on nuclei correspond to small relic abundance, close to our 
lower bound VLh? > 0.025. 

The neutralino cross section on ^''Al was compared with a rough estimate of the 
expected sensitivity of the first stage of the CRESST which uses a low energy threshold 
sapphire cryogenic detector. We have shown that the expected sensitivity to low mass dark 
matter particles of the CRESST experiment can be useful in probing the light neutralino 
scenario (M^ > 2 GeV) predicted by the present analysis based on the unconstrained 
minimal supersymmetric standard model. 
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Figure captions 



Fig. 1: Scanned parameter space a) without and b) with grand unification constraints. 
The existing hmits from accelerator experiments are incorporated. The contour fine 
corresponds to the lower limit on the chargino mass M^± = 45 GeV evaluated for 
tan f3=l (solid lines) and tan (3=50 (dotted lines). The horizontal line at M2=50 GeV 
in the GUT case corresponds to the experimental bound from the gluino mass. The 
upper left and right corners in the GUT case are empty because in these regions 
> 100 GeV. 

Fig. 2: Relic abundance versus neutralino mass in the parameter region allowed by the 
experimental constraints for the four sets B-E listed in Tab. ^ The plot for set A is 
very similar to that for B. The cosmologically interesting relic abundance 0.025 < 
Qh"^ < 1 is within the two horizontal lines. In the case of sets A-C, the upper bound 
in the neutralino mass is given by the slepton mass since we require < Msi- 

Fig. 3: As in Fig. |]with the grand unification constraints. 

Fig. 4: Gaugino purity Zl-^ + d) and bino purity b) for neutralinos in the cosmo- 
logically interesting region. The data sets A-E are combined. The points at bino 
purity Zf-i^ = cos^ 9w = 0.77 correspond to pure photinos. 

Fig. 5: As in Fig. ^with the grand unification constraints. 

Fig. 6: Higgsino composition for the cosmologically interesting neutralinos without grand 
unification constraints. The data sets A-E are combined. The points with Zf^ + 
Zf^ = 1 correspond to pure higgsinos. 

Fig. 7: Relic abundance versus bino purity Zf^ for data sets B and D without grand 
unification constraints. The cosmologically interesting relic abundance 0.025 < 
Qh"^ < 1 is within the two horizontal lines. 

Fig. 8: Bino purity Zf^ versus tan/5 for the cosmologically relevant neutralinos of data 
sets B and D without grand unification constraints. 

Fig. 9: SI versus SD cross sections on ^''Al for the cosmologically relevant neutralinos of 
data sets A, B, C and D without grand unification constraints. 

Fig. 10: SI cross section on ^^Al versus bino purity Zf^ for the cosmologically relevant 
neutralinos of data sets B and D with tan (3 <2 without grand unification constraints. 
The points near Zf^=cos'^ 9 ^ 0.77 correspond to neutralinos with a large photino 
component. 

Fig. 11: SI cross section on ^''Al versus tan/5 for the cosmologically relevant neutralinos 
of data sets B and D without grand unification constraints. 

Fig. 12: SD a) and SI b) cross sections on ^''Al for the cosmologically relevant neutralinos 
obtained without grand unification constraints. The data sets A-E are combined. 
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Fig. 13: SD cross section on ^^Al for the cosmologically relevant photinos and binos with 
purity >0.99. The contribution of the different sets A, B and C is shown. For each 
set, the cross section of pure photino is ~2 times larger than for pure bino. 

Fig. 14: SD cross section on ^'^Al versus bino purity for the cosmologically relevant neu- 
tralinos of data sets B and D without grand unification constraints. 

Fig. 15: SD a) and SI b) cross sections on ^''Al for the cosmologically relevant neutralinos 
obtained for the combined data sets A-E with the grand unification constraints. 

Fig. 16: Cross section ctq (SI+SD) on ^^Na and ^^^I for the cosmologically relevant neu- 
tralinos evaluated without grand unification constraints. The sets A-E are com- 
bined. 

Fig. 17: SD cross section of the cosmologically relevant neutralinos on ^^Ge versus bino 
purity for data sets B and D without grand unification constraints. 

Fig. 18: SI cross section versus the SD part on '^^Ge for the cosmologically relevant neu- 
tralinos of data sets B and D without grand unification constraints. 

Fig. 19: SD a) and SI b) cross section on '^'^Ge versus the neutralino mass for the cosmo- 
logically relevant neutralinos of combined data sets A-E without grand unification 
constraints. 



Fig. 20: As in Fig. ^ with grand unification constraints. 



Fig. 21: Cross section (Tq (SI-I-SD) of the cosmologically relevant neutralinos on ^^''Sn for 
the combined sets A-E without grand unification constraints. 

Fig. 22: Cross section (Jq (SI-I-SD) on ^''Al for the cosmologically relevant neutralinos of 
data sets A-E without grand unification constraints. A rough estimate of the sensi- 
tivity of the first phase of the CRESST experiment is also shown for a measurement 
time of 1 kg-year assuming a flat radioactive background of 1 count /kg/keV/day 
and a detector energy threshold of 500 eV with a full width half maximum energy 
resolution AE=200eV. 

Fig. 23: Normalized interaction rate, in units of count/kg/day/pbarn, showing form fac- 
tor and dark matter flux effects. The rate is evaluated with a zero energy threshold. 
Solid line: rate in a sapphire detector due to interactions on ^^Al. Dashed lines: 
rate in a Nal detector due to interactions on ^^Na and on ^^^I respectively. Dotted 
lines: rate in 100% enriched ^^Ge (upper) and ^^'^Sn (lower) detectors. 

Fig. 24: Cross section (Tq (SI-I-SD) normalized to the reduced mass squared for different 
materials and neutralino compositions. Solid triangles: zino-higgsino mixture with 
40% zino composition. Open triangles: pure higgsino. Solid circles: pure bino. 
Open circles: pure photino. The cross sections are calculated for set B with tan (3=2, 
MA=50GeV and a = —1.24. As can be seen in Fig. ^, zino-higgsino mixtures and 
pure higgsinos are cosmologically relevant above ~20 GeV, while binos and photinos 
are cosmologically relevant above ~2 GeV (without grand unification constraints). 
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Tables 



Tabic 1: Ranges used for the free parameters of the unconstrained MSSM. Calculations with 
the grand unification constraints are performed using the same ranges but with M1/M2 = 0.5 
and M2 > 50GeV. 





mill. 


max. 




M2 


1 


1000 


GeV 


MI/M2 





1 






-1000 


1000 


GeV 


tan/3 


1 


50 





Table 2: Masses in GeV of the squark Mgq, slepton Mgi, and CP-odd Higgs used to 
compute the neutralino relic abundance and the cross section on nuclei for our different trial 
sets. 



set 


Msg 


Msi 


Ma 


A 


45 


45 


25-70 


B 


100 


45 


25-70 


C 


150 


90 


25-70 


D 


200 


200 


200 


E 


1000 


1000 


1000 
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Table 3: Axial coupling coefficients ap(„) for pure photino and pure bino evaluated from 
Eq. (|ll|) using different predictions for the quark spin coefficients AQg. The coupling coeffi- 
cients are shown relative to Mw / Msq='i. 





NQM 


EMC 


Global Fit-1 


Global Fit-2 






0.93 


0.78 


0.80 


0.82 




AQ, 


-0.33 


-0.5 


-0.16 


-0.44 




Ag, 





-0.16 


-0.13 


-0.11 


photino 


dp 


-0.175 


-0.127 


-0.135 


-0.141 




an 


0.020 


0.071 


0.060 


0.054 


bino 


ap 


-0.119 


-0.084 


-0.089 


-0.094 




an 


0.008 


0.045 


0.037 


0.033 
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Table 4: Axial- vector current coefficients Sa for pure photino and pure bino (purity=l) inter- 
actions on the proton-odd nucleus ^"^Al and on the neutron-odd nucleus ''^Ge. The coefficients 
Sa are evaluated using different nuclear models and the nucleon quark spin content coefficients 
listed in Tab. ^. The < Sp(^n) > values used for the Shell Model are unquenched in the case of 
Al [^] and quenched (Q=0.833) in the case of Ge |20|. The coupling coefficients are evaluated 
relative to Mw / Msq=l- 





27A1 (J=5/2) 


73Ge (J=9/2) 




OGM 


Shell Model 


OGM 


Shell Model 


< Sp> 


0.25 


0.3430 





0.011 


<Sn> 





0.0296 


0.23 


0.468 


Sa for pure photino 






NQM 


0.0215 


0.0396 


0.0002 


2.9-10-s 


EMC 


0.0113 


0.0193 


0.0026 


0.0058 


Global Fit-1 


0.0127 


0.0221 


0.0018 


0.0036 


Global Fit-2 


0.0139 


0.0245 


0.0015 


0.0026 


Sa for pure bino 






NQM 


0.0099 


0.0184 


3.3-10-^ 


5.5-10-^ 


EMC 


0.0049 


0.0084 


0.0011 


0.0023 


Global Fit-1 


0.0056 


0.0098 


0.0007 


0.0013 


Global Fit-2 


0.0062 


0.0110 


0.0006 


0.0009 
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Table 5: Axial- vector current coefficients Sa for pure photino and pure bino (purity=l) in- 
teractions on the proton-odd nuclei ^^Na and ^^^I and on the neutron-odd nucleus ^^^Sn. The 
coefficients are evaluated using the OGM for the nuclear structure and the nucleon quark spin 
coefficients derived from a fit of the present data on polarized lepton-nucleon scattering (Global 
Fit-2). The coupling coefficients are evaluated relative to Mw/Msq=l. 





(.1=3/2) 


1271 (.1=5/2) 


"-^Sii (.1=1/2) 


< Sp> 


0.156 


0.07 





< Sn> 








0.261 


Sa for pure photino 


0.0065 


0.0011 


0.0048 


Sa for pure bino 


0.0029 


0.0005 


0.0018 
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